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I. INTRODUCTION 

i— i ! 

The goal of this paper is to extend to multidimensions a new method recently introduced by Fokas [1, 2], as well 
■ as to discuss the relation of this method with the classical methods of images and of integral transforms. It will be 
\ shown that the new method provides an appropriate generalization of the classical approaches, 
f*"*- ■ For simplicity, we will limit our discussion to evolution PDEs on the half line. Evolution PDEs on a finite domain 
are discussed in [3, 4]. The new method can also be applied to elliptic PDEs [5], such as the Laplace [6], Helmholtz [7], 
and biharmonic equations [8] . The extension of this method to nonlinear integrable evolution PDEs is discussed in 

In order to help the reader become familiar with the new method, rather than discussing general initial-boundary 
value (IBV) problems, we will concentrate on the following three concrete, physically significant problems: 



o 

-I— > 



O 



X 



1. The diffusion-convection equation, 



qt=q X x+aq x , < x < oo, t>0, (1.1a) 
q(x, 0) = qo(x), < x < oo, (1-lb) 
q(0,t) = g (t), t>0, (1.1c) 



where a is a real constant. 
2. The linearized Korteweg - de Vries equation (with dominant surface tension), 



qt + q x — Qxxx = 0, < x < oo, t > 0, (1.2a) 
q(x,0) = q (x), < x < oo, (1.2b) 

q(0,t) = g (t), q x (0,t) =gi(t), t>0. (1.2c) 
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3. The multidimensional diffusion-convection equation, 

N 

It = Yffii + a Ai)<J. < ^ < oo, i = 1,2, ... ,7V, t > 0, (1.3a) 

g(xi,X2, . . . ,x N ,Q) = qo(xi,x 2 , . ■ . ,x N ), < Xj < oo, j = 1, 2, . . . , TV, (1.3b) 
<?(0, x 2 , • • • ,xjv,£) = flo^C^, ■ • ■ ,x N ,t), < xj < oo, j = 2,...,N, t > 0, 

g(xi,X2, . . . ,a;jv-i,0,t) = g^\x\,x 2 , ■ ■ ■ ,x N -i,t), < xj < oo, j = 1, . . . , N - 1, t > 0, (1.3c) 

where otj, j — l,...,N are real constants. The functions qo, go, gi, {g$^}i have appropriate smoothness and 
they also decay as x and Xj tend to oo. 

The discussion of the physical significance of the above IBV problems, as well as the derivation of their solution is 
presented in sections II-IV. 

The proper transform in (x,t) 

The proper transform of a given IBV problem is specified by the PDE, the domain, and the boundary conditions. 
For simple IBV problems there exists an algorithmic procedure for deriving the associated transform (see, for example, 
[10, 11]). This procedure is based on separating variables and on analyzing one of the resulting eigenvalue equations. 
Thus, for simple IBV problems in (x, t) there exists a proper x-transform and a proper ^-transform. Sometimes these 
transforms can be found by inspection. For example, for the IBV (1.1) with a = the proper x-transform is the sine 
transform, and if < t < oo the proper ^-transform is the Laplace transform. 

For a general evolution equation in {0 < x < oo, t > 0}, the x-transform is more convenient than the ^-transform. 
For example, looking for a solution of the form exp[ifcx — u)(k)t] in Eqs. (1.1a) and (1.2a), we find that u{k) is given 
explicitly by 

uj(k) = k 2 - iak, u(k) =i(k + k 3 ). (1.4) 

On the other hand, looking for solutions of the form exp[— st + A(s)x] we find that A(s) is given only implicitly, by 

A 2 + aA + .s = 0, A 3 -A + s = 0. (1.5) 

The advantage of the x-transform for an evolution equation becomes clear when the domain is the infinite line 
— oo < x < oo, in which case the x- Fourier transform yields an elegant representation for the initial value problem 
of an arbitrary evolution equation. Furthermore, if an x- transform exists, it provides a convenient method for the 
solution of problems defined on the half line < x < oo. However, in general an x-transform does not exist; this is, 
for example, the situation for the IBV problem (1.2). In this case, until recently one had no choice but to attempt to 
use the i-Laplace transform. For example, Eqs. (1.2) yield 

-q X xx + q~x + sq = qo(x), (1.6a) 
5(0, s) =g (s), q x (0, s) = <h(s), (1.6b) 

where q(x,s), go(s), gi(s), denote the Laplace transforms of q(x,t), go{t), gi(t), respectively. However, this approach 
is rather problematic: (a) The problem is posed for finite t, say < t < T, while the Laplace transform involves 
< t < oo. Thus, the correct transform is 

q(x,s) = / dte~ st q(x,t). 



This yields the additional term — q(x, T) cxp(— st) in the r.h.s. of Eq. (1.6a); some authors use "causality" arguments 
to justify replacing T by oo. (b) If T = oo and if go(t), gi{t) decay for large t, then the application of the Laplace 
transform can be justified. However, since the homogeneous version of (1.6a) involves exp(— A(s)x), where A solves 
the cubic equation (1.5b), the investigation of Eqs. (1.6) is cumbersome. 



The method of images 
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A broad class of IBV problems is often approached by the method of images. Suppose that [12] 

q t + u>(id x )q = 0, < x < oo, t > 0, (1.7a) 

q(x, 0) = q (x), < x < oo, (l-7b) 

g(0,t) = fl (t), *>0, (1.7c) 

where w(-) is an even polynomial of its argument. The method of images then yields a solution as follows. Let 
G(x,x',t) be the Greens function that satisfies 

G t + ui(id x )G = 0, — oo < x < oo, t > 0, 

G(x, 0) = £(x — x'), — oo < a; < oo, < x', 



G(0,i) = 0, i>0, 



and let Q(x,t) satisfy 



Q t + uj(id x )Q = 0, — oo < x < oo, t > 0, 

Q(x, 0) = 0, — oo < x < oo, < x', 

Q(0,t) = </o(t), *>0. 

Note that both G(x, x', t) and Q(x, i) can be readily obtained by means of the Fourier x-transform. Then, the solution 
to Eqs. (1.7) is given by 

q(x,t)= dx' {G(x,x',t) - G(x,-x',t)}q (x') + Q(x,t), x > 0. (1.8) 
Jo 

Indeed, G(x, x', t), G(x, —x', t), Q(x, t) satisfy Eq. (1.7a) and so does their linear superposition, Eq. (1.8). The initial 
condition (1.7b) is satisfied by the term involving G(x, x', t), since G(x, —x', t) and Q(x, t) do not contribute to q(x, t) 
at t = 0. Finally, the boundary condition (1.7c) is satisfied by Q(x,t) alone, since the terms involving G cancel out at 
x = 0. For the latter to be true, it is crucial that oj(-) be an even polynomial of its argument. For example, Eqs. (1.2) 
cannot be treated by the method of images, because lo(1) = i(l + I 3 ) contains only odd powers of I. Neither does the 
method of images apply for Eqs. (1.1), unless a = 0, since lo(1) = —ial + I 2 . 

Sometimes it is possible to apply the method of images after using a suitable transformation. For example, such a 
transformation for Eq. (1.1a) is 

u{x,t) =q(x,t)e ax / 2 . (1.9) 

Using this transformation, Eqs. (1.1) become 

a 2 

ut = u xx —u, < x < oo, t > 0, (1.10a) 

u(x,0) = q a (x)e ax/2 , < x < oo, (1.10b) 
u(0,t) =go(t), t>0. (1.10c) 

Eqs. (1.10) can be solved both by an x-sine transform and by the method of images, provided that w(x,0) decays as 
x — > oo; since the r.h.s. of Eq. (1.10b) involves exp(ax/2) it follows that, for a general initial condition qo(x), this is 
the case iff a < 0. 

The method of images may also work for von Neumann boundary conditions. However, for mixed boundary 
conditions, such as q x (x,0) + (3q(x,0) — go(t), the application of the method of images is far from straightforward. 



Multidimensional IBV problems 

It was noted earlier that IBV problems in {0 < x < oo, t > 0} cannot in general be solved by an x-transform. 
This is the case not only for Eqs. (1.2), but apparently also for Eqs. (1.1) with a > 0. In this case, for problems 
in one spatial dimension one may attempt to use the Laplace transform in t. However, even this approach fails for 
multidimensional problems, since in this case one obtains a PDE in the spatial dimensions, for which there does not 
exist a proper transform. 



The new method 
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FIG. 1: The contour dD+ associated with Eq. (1.1a). 



For equations in one spatial dimension the new method constructs q(x,t) as an integral in the complex fc-plane, 
involving an ^-transform of the initial condition and a ^-transform of the boundary conditions. For equations in 
N spatial dimensions the situation is similar, where the integral representation is now constructed in the complex 
(fci, . . . , fc^-planes. 

For example, for the IBV problem (1.1), it will be shown in section II that this integral representation is 



1 r 1 f 

t) = 7T dk e tkx ^ (k)t q (k) + — dk e lkx -^g(k, t), (1.11) 

27T J-oo 27r JdD+ 

where uj(k) — k 2 — iak, the oriented contour dD + is the curve in the complex fc-plane defined by 



fc/ = f + ^/fc| + (f) 2 . k = k R + ik I: (1.12) 
see Fig. 1, and the functions qo(k), g(k, t) are defined in terms of the given initial and boundary conditions as follows: 

q (k) = dxe~ lkx q (x), Imk < 0, (1.13) 
Jo 

g{k,t) = -q (ia-k)-(2ik + a) [ dr e^ T g ( T ). (1.14) 

Jo 

We note that g(k 7 t) involves go evaluated at ia — k. 



The long-time asymptotics 

The representation obtained by the new method is convenient for computing the long-time asymptotics of the 
solution: Suppose that a given evolution PDE is valid for < t < T, where T is a positive constant. It can be shown 
(see section II) that the representation for q(x, t) is equivalent to the representation obtained by replacing g(k, t) with 
g(k, T). In particular, if < t < oo, then g(k, t) can be replaced by g(k) = g(k, oo). For example, the solution of the 
IBV (1.1) with < t < oo is given by (1.11) with g(k,t) replaced by g(k), 

I* oo 

g(k) = -q {ia-k)-(2ik + a)go(k), g (k) = / dte u( - k ^g (t). (1.15) 

Jo 

Thus, the only time-dependence of q{x,t) appears in the form exp[zfca; — uj{k)t] \ hence it is straightforward to obtain 
the long-time asymptotics, using the steepest descent method. Similar considerations are valid for evolution PDEs in 
multidimensions . 



From the complex plane to the real axis 

In case that the given IBV problem can be solved by an rc-transform, the relevant representation can be obtained 
by deforming the integral representation of the solution obtained by the new method, from the complex fc-plane to 
the real axis. Consider for example the integral in the r.h.s. of Eq. (1.11) for the case that a < 0: It can be verified 
that the functions exp(ikx), qo(ia — k), and exp[w(fc)(r — f)] are bounded in the region of the complex fc-plane above 
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the real axis and below the curve dD + . Thus, the integral along dD + can be deformed to an integral along the real 



axis. 



1 f c 



dk( ,ikx-u(k)t 



(la 



(fc) - q {ia - fc) - {2ik + a) f dr e^ 7 g {r) 

Jo 



(1.16) 



It can be shown that this formula is equivalent to the solution of Eqs. (1.9), (1.10) using the sine transform. However, 
if a > this deformation is not possible. Indeed, qo(ia — k) involves exp[— ix(ia — fc)] = exp[ifc^x] exp[— x(fc/ — a)]; 
this term is bounded for kj > a (and in particular is bounded on dD + ), but it is not bounded in the region of the 
complex fc-plane below the curve dD + . 

We note that even in the cases when it is possible to deform the contour to the real axis, the representation in 
the complex fc-plane has certain advantages. For example, the integral involving the sine transform is not uniformly 
convergent at x — (unless q(0, t) = 0). Also, a convenient way to study the long-time behavior of the representation 
involving the sine transform is to transform it to the associated representation in the complex fc-plane. 



II. THE DIFFUSION-CONVECTION EQUATION 



Physical significance 

Eq. (1.1a) arises, for example, from the diffusion-convection equation 

q T = Dq^ + avq i , a = ±1, < £ < oo, r > 0, (2.1) 

where q(£, r) is a probability density function in the spatial and time variables £ and r, D is a diffusion coefficient 
with dimensions of (length) 2 /(time), and v is a convection field with dimensions of (length) /(time). Eq. (1.1a) is 
the normalized form of (2.1), expressed in terms of the dimcnsionless variables t = (D/v 2 )t and x = (D/v)£. The 
dependence upon the single spatial variable x is justified in systems that are homogeneous in the transverse directions 
to x. The choice a = +1 corresponds to convection (a background drift velocity v) directed toward the origin, 
while a = — 1 corresponds to a background drift away from the origin. In the long-time asymptotic limit convection 
dominates diffusion and determines the fate of the system. Suppose for example that the boundary condition is 
q(0,t) — 0, corresponding to an ideal sink at the origin, for the case, of say, of particle flow. Then, if a = +1 the 
particles will flow to the sink at typical speed v, and the probability density is expected to decay to zero exponentially 
with time. If a = — 1, the particles are drifting away from the sink; in this case one expects a depleted zone near the 
origin, which grows linearly with time, and exponential convergence to a finite level of survival. The case of generic 
boundary conditions is harder to predict by such heuristic arguments. 

The method of images 

We have seen that the method of images is applicable in the singular case of a = 0, and even then the solution 
is straightforward only with pure Dirichlet or von Neumann boundary conditions. For a < the method of images 
can be used after a suitable transform, see Eqs. (1.10), and also provided that the boundary conditions are simple 
enough. For a > 0, the method of images fails to provide a solution for generic initial conditions. We now show how 
the new method serves as a natural extension of the methods of images and of classical integral transforms, even as 
these methods fail. For pedagogical reasons, we expose the new method simultaneously for problems (1.1) and (1.2). 
The physical interpretation of (1.2) is discussed in Section III. 



The new method 

Suppose that a linear evolution PDE in one space variable admits the solution ex.p[ikx — ui(k)t\. For well posedness 
we assume that 



Rew(fc) > 0, for fc real. (2.2) 
The starting point of the new method is to rewrite this PDE in the form 

where 

fu(k)-u(iy 



k-l 



I = -id x . (2.4) 
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For Eq. (1.1a), u(k) = k 2 - iak, thus 



. (k 2 -l 2 -ia(k-l)\ 
X = -i f — ± >-\ q - -[i(k + l) + a]q. 

For Eq. (1.2a), u{k) = i(k + k 3 ), thus 

Hence Eqs. (1.1a) and (1.2a) can be written in the form (2.3), where X is given, respectively, by 

X = -q x - (ik + a)q, X = -q xx - ikq x + (1 + k 2 )q. (2.5) 

(a) The Fourier transform representation 

The solution of Eq. (2.3) can be expressed in the form 

i r°° i r°° 

q(x,t) = — alke lkx -^ t q {k) + — dk e lkx ' u( - k ^ g(k,t), (2.6) 
2tt 7_oo 2tt J_ 00 

where 

poo 

qo{k)= dxe- ikx q {x), lmfc<0, (2.7) 
Jo 

g(k,t)= [ dre w(fe)r X(0,T,/c), k e C. (2.8) 
Jo 

Furthermore, the functions qo(k) and g(k,t) satisfy the global relation 

g(k,t)+q (k)^e^ t q(k,t), lmfc<0, (2.9) 

where q(k, t) denotes the a;- Fourier transform of q(x, t). 

Indeed, using (2.3) it is straightforward to compute the time evolution of q(k,t), 



(e u ^*q{k,tij t = j dx (e-^+^qix^)^ 

= -f dx (e~ ikx+w{k)t x) =e^ k)t X(0,t,k). 

Integrating this equation we find Eq. (2.9). Solving Eq. (2.9) for q(k,t) and then using the inverse Fourier transform, 
we find Eq. (2.6). 

(b) An integral representation in the complex k-plane 

The first crucial step of the new method is to replace the second integral on the r.h.s. of Eq. (2.6) by an integral 
along the oriented curve dD + . This curve is the boundary of the domain D + , 

D+ = [k e C, Imfc > 0, Rew(fc) < 0}, (2.10) 

oriented so that D + is on the left of dD + . For example, for Eqs. (1.1a) and (1.2a) the domains D + arc the shaded 
regions in Figs. 2 and 3, respectively. Each of the curves in Fig. 2 is defined by Eq. (1.12), while the curve in Fig. 3 
is defined by 

ki = ^l + 2,k 2 R . 

Indeed, for Eq. (1.1a), 

Re ui(k) = Re {(&r + ik{) 2 — ia{kn + iki)} — k\ — k 2 + aki. 
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FIG. 2: The domain D+ for Eq. (1.1a). 




i 
i 



FIG. 3: The domain D+ for Eq. (1.2a). 
Thus, D + is the domain of the upper half complex fc-plane specified by 

(*/ - f f - (4 + (f ) 2 ) > o. 

For Eq. (1.2a), 

Rew(fc) = Re{i(k R + ifc/) 3 + i(k R = fc/(fc| - 3fc| - 1). 

Thus, 13+ is the domain of the upper half complex fc-plane specified by 

fcf - 3fc| - 1< 0. 

The deformation of the integral from the real axis to the integral along the curve dD + is a direct consequence of 
Cauchy's theorem. Indeed, the term 

e ikx-u(k)t^ ^ = e ikx f dT e -u{m-T)x(Q, r, k) 
Jo 

is analytic and bounded in the domain E + , 

E + = {ke C, Imfc > 0, Reu(k) > 0}. 

Thus, using Jordan's lemma in E+, the integral along the boundary of E + vanishes. 

For example, for Eq. (1.1a) E + is the domain above the real axis and below the curve dD + ; thus the integral 
along the real axis can be deformed to the integral along dD + . For Eq. (1.2a), E+ is the domain above the curve 
{kj = yfl + 3k R }, thus the integral along this curve vanishes; hence dD + is the union of the real axis and the curve 

{kj = y/1 + 3fc|}. 

(c) Analysis of the global relation 

The second crucial step in the new method is the analysis of the global relation (2.9): This yields g(k,t) in terms 
of qo(k) and the ^-transform of the given boundary conditions. 

Before implementing this step, we note that we have already used Eq. (2.9) in the derivation of Eq. (2.6). However, 
while Eq. (2.9) was used earlier only for real fc, in what follows it will be used in the complex fc-plane. 
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Recall that g(k,t) is defined by Eq. (2.8), where for Eq. (1.1a) X is given by (2.5a). Thus 

g{k,t) = -gi(u(k),t) - (ik + a)g (u(k),t), u(k) = k 2 - iak, (2.11) 

where 

gi (w(k),t) = [ dre^ T q x (0,r), g {u(k),t) = f dr e"^ T g (r); 
Jo Jo 

we have used the notation gi(uj(k),t) and go(u>(k),t) to emphasize that <?i and g depend on k only through uj{k). 
Substituting the expression for g(k,t) from (2.11) into the global relation (2.9) we find 

- gi ((j(k),t) = {ik + a)g {uj{k),t) - q (k) + e u{k)t q{k, t), Imfc < 0. (2.12) 

Our task is to compute g(k, t) on the curve dD + ; since g{k, t) is analytic in D+ this is equivalent to computing g(k, t) 
for k e D + . Eq. (2.11) shows that g(k,t) involves go, which is known in terms of the given boundary condition go{t), 
as well as g\, which involves the unknown boundary value q x (0,t). In order to compute g\ using the global relation 
(2.12), we first transform Eq. (2.12) from the lower half complex fc-plane to the domain D + . In this respect, it is 
important to observe that since g\ depends on k only through uj(k), g\ remains invariant by those transformations 
k — ► v(k) which preserve u)(k). For Eq. (1.1a), the equation u)(k) = u)(u(k)) has one non-trivial root (the trivial root 
is u(k) = k), 

v 2 — iav = k 2 — iak, v{k) = ia — k. 

Let _D_ be defined by 

D- = {k e C, Imfc < 0, Rew(fc) < 0}. (2.13) 

For uj(k) = k 2 — iak, the domains D + and D_ are the shaded regions depicted in Fig. 4. The transformations that 
leave u(k) invariant map the domain {D + ® £)_} onto itself. Thus, if k G D + , then ia — k e _D_. Hence, replacing k 
by ia — k in Eq. (2.12), we find an equation that is valid for k in D + : 

-<?i = -ikg - q (ia - k) + e w(fe) *g(ia - k, t), k e D+. (2.14) 

Replacing — g\ in (2.11) by the r.h.s. of the above equation, we find 

g(k, t) = -{2ik + a)g {uj(k),t) - q {ia - k) + e^ (k)t q{ia - k,t). 

The term e^^q^a — k,t) does not contribute to q(x,t). Indeed, this term gives rise to the integral 

/ dke ikx q(ia - k,t), 

JdD + 

for which wc note: exp(ifcx) is bounded and analytic for Imfc > 0; the term q(ia — k, i) involves exp[— i(ia — k)x] = 
exp(ikRx) exp[— x(kj — a)], which is bounded and analytic for ki > a, i.e., in D + . Thus using Jordan's lemma in D + 
it follows that the above integral vanishes, hence the effective part of g(k,t) is given by (1.14). 

(d) The long-time asymptotics 

Suppose that a given evolution PDE is valid for < t < T, where T is a positive constant. Then, as mentioned 
earlier, we can replace g(k, t) by g(k, T). Indeed, the two representations associated with g(k, t) and g(k, T) differ by 

f dke tkx -^ 1 [ dre^ T X{fd,T,k). 
JdD + Jt 

Since T > t, the coefficient (r — t) of u>(k) is non-negative, thus exp[ikx + w(k)(T — t)] is bounded and analytic in D + , 
and Jordan's lemma implies that the above integral vanishes. 

In the case of Eq. (1.1a) the solution is given by Eqs. (1.11)-(1.13), Eq. (1.15). The long-time asymptotics of q(x, t) 
can be analyzed by the method of steepest descent. For x = vt, the critical point occurs at du/dk = iv, or 

k = —(a + v) = ij. 





a > 



a < 



FIG. 4: The domains D+ and D_ associated with Eq. (1.1a). 

Let v > a > 0, then both integrals on the r.h.s. of (1-11) contribute to the long-time asymptotics, since the path of 
integration can be modified to pass through the critical point in each case. The behaviour of the argument of the 
relevant exponentials near the critical point is given by 

ikvt - uj(k)t = -7 2 < - (A; - ij) 2 t, 

hence the leading-order asymptotics to q(vt, t) is 

1 f°° 1 f 1 

^ J-oo 27T J0D+ 2\/7rf 

The numerical value of qo(ij) and g(ij) are given explicitly in terms of the initial and boundary conditions, see 
Eqs. (1.13) and (1.15). 

For a > v > the second integral on the r.h.s. of (1-11) does not contribute, since the critical point lies outside of 
the domain D + . In this case the leading-order asymptotics is 



q(vt,t) 



1 



2v/7rf 



<?o(*7)e 



For a < and v > — a both integrals contribute and the answer is identical with that of the case v > a > 0. 
Finally, for a < and v < —a neither integral contributes and the leading-order asymptotic behavior is zero. The 
latter is expected on physical grounds, since the (rightward) drift sweeps the probability density q(v i, t) past x = vt. 



III. THE LINEARIZED KORTEWEG - DE VRIES EQUATION 



Eq. (1.2a) is the linear limit of the celebrated Korteweg - de Vries equation: 

qt + q x + Mxxx + Gqq x = 0, A = ±l, (3.1) 

for the case of A = — 1; this corresponds to dominant surface tension. Eq. (3.1) is the normalized form of 



- 2 1 3V 1 3 Th 

or 2 V h a£ 2 3 8 pg 



drj 3 fg d .1 



(3.2) 



where r\ is the elevation of the water above the equilibrium level h, T is the surface tension, p is the density of the 
medium, and g is the free-fall acceleration constant. This equation is the small amplitude, long wave limit of the 
equations describing idealized (inviscid) water waves under the assumption of irrotationality. Eq. (3.1) is obtained 
after transforming to the dimensionless variables t — \J -t-t, x = — a"~ 1//2 £, q = r]/2. 
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FIG. 5: The contours dD^P and dD^ associated with Eq. (1.2a). 

Eq. (3.1) usually appears without the q x term. This is because the Korteweg - de Vries equation is usually studied 
on the full line and then the term q x can be eliminated by means of a Galilean transformation. However, for the 
half-line this transformation would change the domain from a quarter-plane to a wedge. 

Laboratory experiments with water waves typically involve Eq. (3.1) with A = 1, q(x, 0) = 0, and q(0, t) a periodic 
function of t. Thus the linearized version of Eq. (3.1) with A = 1 is valid until small amplitude waves reach the 
opposite end of the water tank. Eq. (1.2a) is valid under similar circumstances, in the case of dominant surface 
tension. 

It is interesting to note that while in the case of A = 1 the problem is well posed with only one boundary condition 
at x = 0, in the case of A = — 1 the problem is well posed with two boundary conditions. The case of A = 1 is solved 
in [2] . Here we solve the case of A = — 1 . 

The solution of the IBV problem (1.2) is given by Eq. (1.11), where ui{k) = i(fc + fc 3 ), qo(k) is the Fourier transform 
of qo(x) (see Eq. (1.13)), g(k,t) is defined by 

g(k, t) = -q {v) + i{v - k) f dr e^ 9l (r) + (k 2 - v 2 ) f dr e^g (r), (3.3) 

Jo Jo 

-k-iV3k 2 + 4: (1) -fc W3fc 2 + 4 (2) 

v = for k £ D\ ' v = for k G D\ ' , 

2 + 2 + 

and dD + is the union of the boundaries of and depicted in Fig. 5 (the relevant curve is {ki = y/Sk^ + 1, kj > 
0}). 

Indeed, the representation (1.11) was derived in section II. Using the definition of g(k,t), Eq. (2.8), and recalling 
that X is given by Eq. (2.5b), we find 

g{k, t) = -ga(u{k), t) - ifc&(w(fc), t) + (1 + k 2 )g (u(k), t), u(k) = i(k + fc 3 ), (3.4) 

where g±, go are the first, second integrals appearing in the r.h.s. of Eq. (3.3) and gi involves the unknown boundary 
value q xx (0,t), 

fc(cj(k),t) = f dTe^ T q xx (0,T). 
Jo 

Substituting the expression for g(k, t) from Eq. (3.4) into the global relation (2.9), we find 

-&(w(fc), t) - t) - (1 + k 2 )g (iu(k), t) - q Q {k) + e^ l q{k, t), Im k < 0. (3.5) 

The equation u)(k) — oj{v(k)) has two nontrivial roots, 

+ = fc + fc 3 , v 2 + kv + k 2 + l = Q, v h2= ~ kTW 3 t-±± ; (3.6) 

i^(jfe) - e 4 "/ 3 /c, v 2 {k) - e 2m ' z k. s k^oo. 

The domains D + and D_ (see Eqs. (2.10) and (2.13)) are the shaded regions depicted in Fig. 6a, where the two 
relevant curves are kj = ±y3H + l; the limit of the domains as k — > oo are depicted in Fig. 6b, where each of the 
relevant angles is 7r/3. 

If ^ < arg k < 4f then ^ < arg k + ^ < 2tt, and 2tt < arg k + 4f < ^. Thus 

k G Df , v 2 (fc) G D_ , vy (k) G DP . 
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(a) 




FIG. 6: The domains D+ and D_ for Eq. (1.2a). 



Similarly, 

k e D ( +\ u 2 {k) e v x {k) e r>_. 

Thus, replacing fc by v\{k) and 1*2 (fc) m (3-5), we find 

-g 2 - lv x g x - (1 + ^ 2 ).9 - gb(i^) + e w ( fc ) t g(i/ 1 (fc), i), k G D* 1 '; 

-.92 = - (1 + i/I)flO - <7o(^ 2 ) + e"( fc )*q(^(fc),<), fc G D^. 
Replacing —§2 in (3.4) by the r.h.s. of the above equations, we find 

g{k, t) = i{u - k) 9l + {k 2 - u 2 )g Q - q Q {u) + e^qiy, t), k G D+. 

Due to analiticity considerations the term q(v,t) does not contribute to q(x, t), see section II; thus the effective part 
of g(k,t) is given by (3.3). 

The long-time asymptotics 

The long-time asymptotics can again be analyzed by replacing g(k,t) with g(k } 00) = g(k) and following standard 
steepest descent or stationary phase expansions. One thus finds [13]: For v = x/t > 1, q(vt,t) satisfies 

q(vt, t) = [(P( 7 ) - P{v(rt))e* + (P(- 7 ) - P{v{-i))e^] + 0(t" 3 / a ), 

as t — > 00, where 
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v-l 



3 ' 

<p(t) = 2 7 3 t-7r/4, 
v(k) = -i(fc + iV3fc 2 + 4), 
P(fc) = <7o(fc) + (fc 2 + l) 3o (fc)-^. 9l (fc), 

qo(fc) is defined in terms of the initial condition (Eq. (1.15)), and gi(k), go(k) are the first, second integrals appearing 
in the r.h.s. of Eq. (3.3) evaluated with t — > 00. For v < 1, q(vt,t) decays faster than any algebraic power of i, as 
i — > 00. 
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IV. THE DIFFUSION EQUATION IN MULTIDIMENSIONS 



Physical significance 

This problem is the generalization of problem (1.1) to multidimensions. It arises from the diffusion-convection 
equation 

JV 

'/ X 1 ^'/; i ' O'/; '• < & < oo, r > 0, j = l,2,...N, (4.1) 

which represent isotropic diffusion in a quadrant of TV-dimensional space, with diffusion coefficient D, and under a 
background convection (drift) field (vi,V2, ■ ■ ■ , vn). Eq. (1.3a) is obtained upon passing to the dimensionless variables 
t = (D/v 2 )t and Xj — (D/v)£j, j = 1, 2, . . . , TV, where v is a typical speed such as the rms v = (J2j vf) 1 / 2 . 

The new method 

For the solution of the IBV problem (1.3) we introduce the following notations: 

N N 

• x = (x 1 ,...,x N ), k = (fei, . . . , k N ), k-x = ^kix u fc (j) • x {j) = ^ k t x t . 

l l=U& 



p pOO pC 

/ dk = / dh--- 

JR N J-oo J-< 



dk 



N- 



Cj is the curve {(%) 7 > 0, (fy)/ = ^ + ^[(fc,)«] 2 + (^) 2 , j = 1, 2, ... , TV}. 
v j =ia j -k j , j = 1,2, ...,7V. 

AT 



pOO pOC 

• Qo{k)= dxi--- dx N e~ ik ' x q {x). 

Jo Jo 

rt />oo roo j-oo ^oc 

• g j) (k,t) = / dr / dx 1 --- dxj-i dx j+1 --- dx N e~ ik3 ' x 3 +uj(k)T g { ^ {x-l, . . . ,x j -- L ,x j+1 , . . . ,x n ,t), 

Jo Jo Jo Jo Jo 

i = 1,2,..., TV. 

AT 

• Q(fc, t) = «o(fc) - + ctj)$\k, t). 

i 

The solution of the IBV problem (1.3) is given by 
g(x,t)= / dfce*— W*,?o(fc) 

( 27r ) JR» 

1 ^ pOO pOO p pOO pOO 

-T^nY, dk!--- dkj-i dk, dk j+1 --- dk N e ik - x - u ^ t (ik j +a j )g^ ) {k,t) 

-J^N (yf c dk l dfc 2 ' ' ' ^e^"^'^!, fc 2 , . . . , fcjy) + Cp) 

+ (2^p (X dh Ic dk2 I ^ " ' I dkNelk ' X ~ U(k)t ®^ ■■■' k *) + CP^j 

f dh f dk 2 --- f rffcjve^-^Q^i,^,...,^), (4.2) 



_... + (zil 

(27T) 



where CP denotes cyclic permutations (i.e., all possible relevant combinations). 
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For the derivation of Eq. (4.2) we introduce the additional notations: 

• = (x\, . . . ,Xj-i,x j+ i, . . . ,x N ), k^ = {ki,...,kj-i,k j+ i,...,k N ). 

f- f-OC f-OC 

• / dx = dxi - I dxN- 
Vr« Jo Jo 

f- f-OC f-OC f-OC f-oc 

• / dx^ = / dxi ■ ■ ■ dxj-i / dxj + i ■ ■ ■ dxjy. 

JR^' 1 JO JO Jo Jo 

f- f-oc f-oc f-oc f-oo 

• / dk^ = / dki ■ ■ ■ dkj^i / dkj + i ■ ■ ■ dkN- 

J~R(.3) J — oo •/ — oc J— oo J —oc 

• e = e ik ' x ~ u '( k ) t 

• E = e w(k)t . 

Substituting exp[ifc • x — w(fc)i] in Eq. (1.3a) we find that ui(k) is given by the expression defined in the notations. 
Eq. (1.3a) can be written in the form 



N 

^-ik-x+^t^ ^ + ^ (e-^+^Xj^ = 0, (4.3) 



i 

where 

Xj = -q Xj - (ikj + aj)q, j = 1,2, . . . , TV. 
Let q(k,t) denote the TV-dimensional Fourier transform of q(x,t). Using Eq. (4.3) we find 

(e^Wqik, tfj t = J dx (e- ik - x+ ^ k)t q(x, i)) ^ 

= -Y dxA e -*'+»%. = V e-^+^Xjl^o. 

~,Jo V ' x i ~~f 



J = 1 J0 ~ J 3=1 



Integrating this equation with respect to t, we obtain 

N 



q (k) + J29 {i) (k,t)=e^ t q(k,t), (4.4) 
l 

where qo{k) is the TV-dimensional Fourier transform of qo(x) and g^\k, t) is defined as follows: 

g^(k,t) = - fdr f dx^ e- ikU) - xU)+ ^ k >[q Xj (x^ ,t) + (ikj + aj )q(x^ ,r)\, 
Jo Jr? 



or 



g^(k,t) = g[ j) (-ik^,oj(k),t) - (ikj+a^i-ik^Mk^t), j = 1,2, . . . ,N, (4.5) 

where g^ is defined in the notations and g[^ denotes the integral of — q Xj . 

Solving Eq. (4.4) for q(k, t) and using the inverse TV-dimensional Fourier transform, we find 

N 

"■-»<*>'[9<»-(i* j+<tj )5»]. (4.6) 



t) = WfL" k e """" < '"'*w * « 



The terms and contain fcj only through cj(fc), thus the relevant integral in the complex fc^-plane can be 
deformed from the real axis to the curve Cj (see the discussion in section II). Hence Eq. (4.6) yields the first two 
integrals appearing in (4.2) plus the term 

<^%L dkU) L^ e "^- (47) 
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Using the global relation (4.4), it can be shown that the above term yields the remaining expressions appearing in 
(4.2). For pedagogical reasons, we give the details for N = 3: In this case the global condition becomes 

g[ 1 \-ik 2 , -ik 3 ) + gf \-ik u -ik 3 ) + g[ 3 \-ik 2 , -ik 3 ) = -Q{h,k 2 , k 3 ) + Eq(k 1 ,k 2 ,k 3 ), lmkj < 0, (4.8) 

where E, Q are defined in the notations and, for simplicity of notation, we have dropped the t- and w(/c)-dependence. 
Since g[^ in Eq. (4.7) is integrated along C\, we replace k\ by v\ in Eq. (4.8) and then solve the resulting equation for 
g^\—ik 2 ,—ik 3 ). Similarly, for g^\ g^\ we replace k 2 , k 3 in Eq. (4.8) by v 2 , v 3 . Thus the expression (4.7) becomes 

- / dk 2 dk 3 dki e[gf\-iui,-ik 3 ) + g^\-iui, -ik 2 )] 
Jr 2 Jd 

- I dkidk 3 dk 2 e[g^\-iv 2 , -ik 3 ) + g^\-iki, -iv 2 )} 
Jr2 Jc 2 

-[ dk x dk 2 ( dk 3 e[g ( ^{-ik 2 ,-iu 3 )+gf\-ik u -iu 3 )] (4.9) 

JR 2 Jc 3 

- / dk 2 dk 3 / dk\ eQ{v\,k 2 ,k 3 ) - I dk x dk 3 I dk 2 eQ(ki,v 2 ,k 3 ) - / dk x dk 2 I dk 3 eQ{ki,k 2 ,v 3 ). 

JR 2 Jd JR 2 JC-2 JR 2 JC 3 

In addition to these terms we obtain three terms involving q, but these terms vanish due to analyticity considerations; 
for example, one of these terms is 

/ dk 2 dk 3 dkie ihx q(i'i,k 2 ,k 3 ,t), 
Jr 2 Jd 

which vanishes since both e lk ' x and q(fi,k 2 ,k 3 ,t) are bounded and analytic in the region of the complex fci-planc 
above the curve C\. 

The last 3 terms in (4.9) are part of the expression for q(x,t); regarding the first 3 terms in (4.9), we note the 
following: the two integrals involving g^ contain k 2 only through u>(k), thus for these integrals the contour with 
respect to k 2 can be deformed from the real k 2 -&xis to the curve C 2 ; similarly for the integrals involving and g^ . 
Thus we obtain 



/ dk 3 dk 2 / dki e[g^\-w 2l -ik 3 ) + gf\-ivi,-ik 3 )], 
Jr Jc 2 Jd 



(4.10) 



plus two more similar integrals, which can be obtained from (4.10) by cyclic permutation 1 — > 2 — » 3 — » 1. The square 
bracket appearing in (4.10) can be expressed in terms of g^ using the global relation evaluated at ki — > ui, k 2 — ► v 2 , 

-g^\-iv 2 , -ik 3 ) - g^\-iui, -ik 3 ) = g[ 3 \-ivi,-iv 2 ) + Q(vi,v 2 , k 3 ) - Eq(vi, u 2 ,k 3 ). 

Due to analyticity considerations the term involving q vanishes, thus the expression (4.10) together with the other 
two similar expressions yield 

/ dk 3 I dk 2 I dkieg^[\—iv-L,—iv 2 )+ I dk 2 / dk 3 / dk\ eg^\—ivi, —iv 3 ) 

J R J C2 J Ci J R J O3 J C\ 

+ dki dk 2 / dk 3 eg^\-iv 2 , -iv 3 ) 

J R J C*2 C3 

+ / dk 3 / dk 2 / dk\eQ{v\,v 2 ,k 3 ) + / dk 2 / dk 3 / dk x eQ(v u k 2 ,v 3 ) 

J R J C2 ^ Ci J R J O3 J Ci 

+ / dk x I dk 2 I dk 3 eQ(ki,v 2 ,v 3 ). (4.11) 

J R J C2 •* C3 

The last three terms above, are part of the expression for q(x, t); regarding the first three terms we note the following: 
the integral involving g^ contains k 3 only through w(fc), thus the integration along the real k 3 -axis can be deformed 
to the curve C 3 ; similarly for the other two integrals. Thus the first three terms of (4.11) yield 

/ dki dk 2 dk 3 e\g^\-iv 2 ,-iv 3 ) + g^\-iv 1 ,-iv 3 ) + g[ 3 \-ivi,-iv 2 )]. (4.12) 
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Using the global relation evaluated at k\ — > u\, k 2 — > z/2, &3 — ► ^3, the square bracket appearing in (4.12) can be 
replaced by 

-Q{V\,V2,V 3 ) + £g(l/i,I/ 2 ,I/3, *)• 

Thus the expression (4.12) yields 

J C\ J Co. J O3 

plus a term involving q that vanishes due to analiticity considerations. 

The long-time asymptotics 

The long-time asymptotics may be analyzed by the method of steepest descent, after replacing the variables (k, t) 
with gW(fc) =gW(k,oo). 
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